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Abstract 

In this paper another class of Dark Matter candidate particles: the pseu- 
doscalar and scalar light bosonic candidates, is discussed. Particular care is 
devoted to the study of the processes for their detection (which only involves 
electrons and photons/X-rays) in a suitable underground experimental set-up. 
For this purpose the needed calculations are developed and various related as- 
pects and phenomenologies are discussed. In particular, it is shown that - in 
addition to the WIMP cases already discussed elsewhere - there is also possi- 
bility for a bosonic candidate to account for the 6.3 a C.L. model independent 
evidence for the presence of a particle DM component in the galactic halo ob- 
served by DAMA/Nal. Allowed regions in these scenarios are presented also 
paying particular care on the cosmological interest of the bosonic candidate. 

Keywords: Dark Matter; axion-like particles; light bosons; underground Physics 
PACS numbers: 95.35. +d, 14.80.Mz, 29.40.Mc 



1 Introduction 

In order to investigate in a model independent way the presence of a DM particle 
component in the galactic halo we have exploited the effect of the Earth revolution 
around the Sun on the DM particles interactions on suitable underground detectors. 
In fact, as a consequence of its annual revolution, the Earth should be crossed by a 
larger flux of DM particles in June (when its rotational velocity is summed to the one 
of the solar system with respect to the Galaxy) and by a smaller one in December 
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(when the two velocities are subtracted). This offers an efficient model independent 
signature, able to test a large interval of cross sections and of halo densities; it is 
named annual modulation signature and was originally suggested in the middle of 
'80 by [1]. The annual modulation signature is very distinctive since a DM particle 
induced seasonal effect must simultaneously satisfy all the following requirements: the 
rate must contain a component modulated according to a cosine function (1) with one 
year period (2) and a phase that peaks roughly around ~ 2 nd June (3); this modulation 
must only be found in a well-defined low energy range, where DM particles induced 
events can be present (4); it must apply to those events in which just one detector 
of many actually "fires", since the dark matter particle multi-scattering probability 
is negligible (5); the modulation amplitude in the region of maximal sensitivity must 
be <7% for usually adopted halo distributions (6), but it can be larger in case of 
some possible scenarios such as e.g. those in refs. [2, 3, 4]. Only systematic effects 
able to fulfil these 6 requirements and to account for the whole observed modulation 
amplitude could mimic this signature; thus, no other effect investigated so far in the 
field of rare processes offers a so stringent and unambiguous signature. 

The DAMA/Nal set-up has exploited such a DM annual modulation signature 
over seven annual cycles [5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. It has achieved a 6.3 a 
C.L. model independent evidence for the presence of a DM particle component in the 
galactic halo. Some of the many possible corollary quests for the candidate particle 
have been carried out mainly focusing so far the class of DM candidate particles 
named WIMPs [6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. In literature several candidates for 
WIMPs have been considered - all foreseen in theories beyond the Standard Model 
of particle Physics - such as the neutralino in the supersymmctric theories [16, 17, 
18], a subdominant sneutrino (the spin-0 supersymmetric partner of the neutrino), 
a subdominant neutrino of the fourth family (see the analysis including the results 
of the indirect search in ref. [19]), a dominant sneutrino in supersymmetric models 
with violation of lepton number, where two mass states and a small energy splitting is 
present, as reported in ref. [2], the particles from multi-dimensional Kaluza-Klcin-like 
theories, etc. In addition, other possibilities still exist with a phenomenology similar 
as for the WIMP cases: the mirror Dark Matter particles [20] , the self- interacting dark 
matter particles [21], etc. Moreover, in principle even whatever particle with suitable 
characteristics, not yet foreseen by theories, can be a good candidate as DM in the 
galactic halo. 

In the present paper the class of light bosonic candidates - either with the pseu- 
doscalar or with the scalar coupling will be investigated in some of the possible 
scenarios. Moreover, also here we will not yet consider either the case of halos with 
caustics or with streams contributions [4, 3], which can have significant impact for 
all the DM candidate particles; we plan to exploit these scenarios in near future by 
further devoted analyses. 

It is worth to note that the direct detection process for light bosonic DM candidates 
is based on the total conversion in Nal(Tl) crystal of the mass of the absorbed bosonic 
particle into electromagnetic radiation. Thus, in these processes the target nuclei recoil 
is negligible and is not involved in the detection process; therefore, signals from these 
light bosonic DM candidates are lost in experiments based on rejection procedures of 
the electromagnetic contribution to the counting rate. 
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2 The light bosons candidates and their main inter- 
actions with ordinary matter 

Let us firstly remind the Peccei-Quinn (PQ) axion; it is a pseudoscalar particle hypoth- 
esized in order to solve the strong CP problem. Several models exist which describe the 
coupling of such a particle to ordinary matter; the most popular ones are the DFSZ, 
where the particle is directly coupled with electron [22], and the KSVZ, where such a 
coupling is not present at tree-level [23]. However, many models exist for axion-likc 
particles, that is particles having similar phenomenology with ordinary matter as the 
axion, but which allow values for the coupling constants and for the mass significantly 
different from those foreseen in the DFSZ and KSVZ models. For example, we mention 
the axion itself in the Kaluza-Klcin theories [24] , where it would have similar couplings 
as in the DFSZ and KSVZ models, but much higher mass states or the "exotic" ax- 
ion models proposed by [25]. Other candidates are pseudo-Nambu-Goldstone bosons 
related to spontaneous global symmetry breaking different from the U(l)pg hypoth- 
esized by Peccei-Quinn, such as the pseudoscalar familon in the case of the family 
symmetry or the Majoron for the lepton number symmetry [26, 27]. 

The previous cases mainly refer to light bosonic candidates with pseudoscalar cou- 
pling; obviously, cases with scalar coupling can be considered as well, such as e.g. the 
scalar familon [28] and the sgoldstino [29]. 

Moreover, it is also interesting to remark that some indirect astrophysical obser- 
vations: i) the Solar corona problem; ii) the X-rays flux detected by ROSAT in the 
direction of the dark side of the Moon; iii) the X-rays background radiation in the 
2-8 keV region measured by CHANDRA (XRB); iv) the excess of X-rays from clus- 
ters of galaxies; have recently been analysed in a model of axion-like particles with 
mass in the keV range and coupling to photons g ail of the order of 10 _13 GeV _1 [30], 
that is requiring a model with expectations for the coupling constants and masses well 
different with the respect to those expected in the DFSZ and KSVZ models. 

It has also been argued that the existence of axion-like particles may account for 
the high energy cosmic rays [31]. Also in this case the bosonic candidate particle can 
have masses of few keV and g ail of the order of 10~ 15 GeV -1 , that is still couplings 
well different than those expected by the DFSZ and KSVZ models. 

Finally, a keV Majoron has been suggested as DM particle [32] and a ~ keV DM 
pseudoscalar candidate has also been taken into account in ref. [25, 33]. 

In this paper, we consider a light bosonic candidate, either with pseudoscalar or 
with scalar couplings, of ~ keV mass as DM component in the galactic halo 2 . Several 
mechanisms can be advocated for the production of these particles in the early Universe 
(see e.g. ref. [32, 25, 33]); in particular, it has been demonstrated that these particles 
can be of cosmological interest, providing many configurations with il <~ 0.3. Moreover, 
such a DM keV candidate can also be able to explain the galactic scale in the structure 
formation problem [32] . 

The diagrams of Fig. 1 show the main processes involved in the detection of a DM 
light bosonic particle (here generically named a) both in the pseudoscalar and in the 

2 In fact, considering the phenomenology of this candidate (see later), a keV-scale bosonic candidate 
naturally arises as an additional solution to the observed model-independent DAMA/Nal annual 
modulation signal. 
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scalar interaction types. 




Z,e" 



2,e 
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or Photoelectric-like 



Primakoff 



Figure 1: Possible diagrams for the direct detection processes of a light boson particle. 
On the left: "Compton - like" effect; on the center: "Axioelectric" or photoelectric-like 
effect; on the right: Primakoff effect. See text. 

As first, let us recall that in the case of interest here the DM light bosonic candidates 
are non-relativistic since they should be trapped in the galactic halo 3 . 

On the left of Fig. f the "Compton - like" effect is shown; in fact, considering that 
the light boson particle couples to charged fermions of the ordinary matter (either 
electrons or u and d quarks), it is absorbed as in the ordinary Compton effect and it 
is emitted a photon, whose energy is ~ m a if - as in the case of interest here - m a is 
much lower that the charged fermions masses. 

In the center of Fig. 1 the "axioelectric effect" or photoelectric-like effect is shown; 
there the a particle - in the hypothesis of a non-zero coupling to electrons - is absorbed 
by a bound electron of the atom exciting or ionizing it. Thus, the sum of the kinetic 
energy of the photoelectron and of the energy produced by X-rays and Auger electrons 
in the re- arrangement of the atomic shell is ~ m a ; this is the detected quantity. This 
effect is phenomenological similar to the ordinary photoelectric effect. 

On the right of Fig. I the Primakoff effect is shown; in this case a is coupled 
with two photons through a charged fermions loop (leptons or quarks). One of the 
two photons is virtual and is related to the electric field of the atoms in the Nal(Tl) 
crystal, the other photon has energy ~ m a and can be detected. 

For the sake of completeness, we note that the elastic scatterings of these light 
bosonic particles on atoms, nuclei and electrons are not considered here since their 
effect is negligible in underground direct detection experiments. 

In all the processes described above the total (including the secondary processes: 
X-rays and Auger electrons) energy release, E re i, in the detector (providing that its 
detection efficiency is ~ 1 for low-energy electrons and low-energy photons) matches 
the total energy of the a particle, E a ~ m a since the a velocity is of the order of I0 _3 c. 

3 Moreover, in case these particles are relics from the early Universe, they should have been cooled 
down by the cosmological redshift; thus, they can be non-relativistic now, similarly as the case of 
massive neutrinos. Other scenarios for their production and for their non-relativistic properties can be 
in principle considered as well. On the other hand, in order to induce the observed annual modulation 
signal, they cannot obviously be relativistic (see also later). 
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Therefore, the differential cross section can generally be written as: 

= cT tot (v)5(E re i - E a ) ~ a tot (v)5(E re i - m a ) (1) 

aHirel 

where a to t generally depends on the DM particle velocity in the laboratory frame, 
\v\ = v. 

Taking into account the energy resolution of the detectors, A, the total differential 
counting rate per unit target as function of the detected energy, E ee , is: 



(E„ -moVf 



d£ e e - J "^e/ V ^: A J a V d 3 v dErel - tltot X e 2A (2 j 

where i? tot is the area of the gaussian peak centered at the a particle mass. 

The incoming flux of a particles is given by d(j>(v) — ^f{v + v m )d 3 v, with p a = 
£ ■ Phaio local density (in GeV • cm~ 3 ) of the a particles in the galactic halo (phaio is 
the local halo density and £ < 1 is the a particle fraction amount of the local density). 
In the following for simplicity we will assume £ = 1, since the scaling of the presented 
results can be derived straightforward for the other cases. Moreover, f(v + v®) is 
the DM particle velocity distribution, which also depends on the Earth velocity in the 
galactic frame, v®, and, thus, on the time along the year. Therefore, the total counting 
rate per unit target is: 

Rtot = [ cPv^fiv + v^atotiv) (3) 
J m a 

and, generally, it can depend on the time along the year. As an example, we cite here 
a case, that we will see again in the following: let us assume that a to t(v) oc v; thus: 

R tot cx J d 3 v v 2 f{v + v m ) = (v 2 ) (4) 

Defining the velocity of a particle in the galactic frame as v g = v + , we obtain for 
non rotating halo: (v 2 ) — (v 2 ) + v®. Since the Earth velocity in the galactic frame, , 
is given by the sum of the Sun velocity, Vq , and of the Earth's orbital time-dependent 
velocity around the Sun, vsE{t), we obtain neglecting the v 2 SE term: 

(v 2 ) ~ (Vg) + Vq + 2v G ■ vsE(t) ^ (v 2 ) + Vq + VQV S ECos(u(t - t )) (5) 

In the last equation we have considered the angle of <~ 60° of the terrestrial orbit 
with the respect to the galactic plane, u = 2n/T with T = 1 year and the phase to 
corresponding to ~ 2nd June (when the Earth's speed in the galactic frame is at the 
maximum). The Sun velocity can be written as |tT©| ~ vo + 12 km/s, where vq is the 
local velocity, whose value is in the range 170-270 km/s [8, 34]; the Earth's orbital 
velocity is vse — 30 km/s. Moreover, (v 2 ) depends on the halo model and on the Vq 
value (just for reference, in the particular unphysical case of isothermal halo model: 
(v 2 ) = |i>o ). Finally, in this illustrative case R to t = So + S m ■ cos(uj(t — to)), giving a 
time dependent contribution with amplitude depending on the adopted halo model. 

In the following, the contributions of each process and of each coupling to the 
&tot(v) will be shown, as well as the contributions to R to t through the cq. (3). Note 
that in all the following formulae we assume % = c = 1. For clariteness, hereafter 
the DM bosonic candidate particle will be named either a if it has a pseudoscalar 
interaction or h if it has a scalar one. 



5 



3 The pseudoscalar case 



In Fig. 2 the diagrams, which describe the interactions of the pseudoscalar a particle 
with the fermions / (in general leptons and quarks) and the "anomalous" coupling to 
two photons, are shown. 




Figure 2: Coupling of the pseudoscalar a particle to the fermions / (a case) and to the 
photons through the loop of charged fermions (b case) . 



The effective Lagrangian of interaction of the pseudoscalar a particle is: 



Lint = i g a f f a /7 5 / + 



4 aF^F*" 



(6) 



where the second term describes the effective coupling to the photons through a charged 
fermion loop. As a consequence the pseudoscalar a particle can decay into fermions 
and photons. The decay in couples // could obviously be possible if a would have a 
mass larger than two times the fermion mass, m/. 

In order to be of cosmological interest as DM candidates in the Universe the lifetime 
(t ) of a should roughly be of the order or larger than the age of the Universe 4 , although 
this condition might be released in case other exotic scenarios would be considered. 
Its lifetime can be derived considering just the second term of cq. (6) under the 
assumptions that m a < 2m e ~ 1 MeV and that the possible decay into massive 
neutrinos (whose coupling constants g ai)v can be expected to be very low) would be 
ncglible with the respect to the two-photon decay contribution: 



64tt 



fJa 



(7) 



77 



The effective coupling constant to photons g ail can be calculated by taking into 
account the charged fermions loop with Qf charge. For a pseudoscalar particle having 
m a « rrif one gets [35]: 



E 



9affQ) a 

irrrif 



m a I ml 



(8) 



4 The depletion of these particles over the course of the history of the Universe is largely dominated 
by the decay processes. Therefore, the effect of depletion through conversion processes on ordinary 
matter are consequently not addressed here. 
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3.1 Direct detection of the pseudoscalar particle 
3.1.1 The Compton-like effect 

Let us take into account the diagram of Fig. 1, where a pseudoscalar particle having 
a momentum k = m a v is absorbed by a charged free and point-like fermion and a 
photon is emitted. In this case, in the non-relativistic limit, one gets (see Appendix): 



>»:, 



considering that the fermion recoil is negligible. 

Since the pseudoscalar interaction gives a spin-dependent contribution (see Ap- 
pendix), the coherent conversion on closed (spinless) shells vanishes. Hence, the 
counting rate Rc-i expected in a detector, made of various species of target fermions 
(electrons and nuclei in Nal(Tl)) identified with the / index, is: 

Rc-i = N T ]T ±^J— Ff { q ) (10) 

/ ^ 

where N T is the number of target atoms (either Na or I atoms in Nal(Tl) detectors) 
and the sum runs over all the electrons and over the nuclei (if they have unpaired 
nucleon, as it is the case for 23 Na and 127 I nuclei). Ff is the incoherent scattering 
function [36] of the transferred momentum, q ~ m a . 

As far as regards electrons, the incoherent scattering function takes into account 
their binding in the atoms. In particular, a good approximation for non relativistic a 

2 

particle is to consider that Ff(m a ) = 1 when the electron recoil energy (~ 2 (m"+m ) ) ^ s 
larger then the binding energy of the bound electron, zero elsewhere. For this purpose, 
we report in Tab. 1 the binding energies for the electrons in Na and I atoms. 

For the case of nuclei, eq.(10) still holds when using for Qf the effective charges of 
the nuclei, Z e ff, (which take into account the screening of the electrons) and for the 
coupling constants the following scaling law (valid for nuclei with unpaired proton, as 
Na and I are): 

9aAA _ 9app ^ ^ / ^ 9auu 9 add j (H) 



m,A rn p 3 \ m* m_ 

where m* and m* d are the constituent up and down quark masses. This scaling law has 
been determined with a similar procedure than those used for the evaluation of nucleon 
magnetic moments in NRQM [38]. The incoherent scattering function for the nuclei 
can be considered as one for the energy of interest. As it can be seen, the contribution 
of the nuclei to the total cross section can be generally neglected because of the larger 
mass of the nuclei with the respect to the electron mass. 

Finally, it is worth to note that, due to the dependence of crc-l,f on the inverse 
of the a velocity, in this particular case Rc-i does not depend on the time along the 
year (S m = 0). Therefore, this process only contributes to the constant part, So, of 
the signal. 
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Table 1: Binding energies of elettrons in 23 Na and 127 1 [37] 



Shell 


E Na (eV) 


Ei (eV) 


Is 


1070.8 


33169 


2s 


63.5 


5188 


2P1/2 


30.4 


4852 


2^3/2 


30.5 


4557 


3s 


few eV 


1072 


3P1/2 




931 


3j?3/2 




875 


3rf 3 /2 




630.8 


3^5/2 




619.3 


4s 




186 


4pi/ 2 




123 


4^3/2 




123 


4d 3 / 2 




50.6 


4rf 5 /2 




48.9 


5s 




few eV 


5p 




few eV 



3.1.2 The axioelectric effect 



The process similar to the usual photoelectric effect has been investigated. Hence, 
it has been calculated (see Appendix) the differential cross section for the process of 
absorption of the pseudoscalar a particle (having a momentum k) by an atomic electron 
(described by the wave function ip n i m ) which is extracted (with final momentum p) 
from the atomic shell with quantum numbers nlm: 



d(JA-e .nlm 2 

dn ~ 9age 



MM \p\ 

2m e (2tt) 3 



-i(p—k)x 



1pnlmd 3 X 



(12) 



Therefore, since in this case q = p — k ~ p, we expect <JA-e,nlm oc v and a non- 
zero modulation term appears, as described before. To evaluate the ip n im Fourier 
transform integral, we have approximated the wave functions to the hydrogenic case: 
/ e-^ nlm d 3 x = (2Tr) 3 / 2 F n [(q) ■ Y; TO (f2), where Y"; m are the spherical harmonics. 
Defining Q = -^L with qo = l/ fl o an d a o Bohr radius, one gets [39]: 



Fni(q) 



2(n-l-l)\ 



1/2 



'^^" pgrp e 1 ,-! (^o)~ 3 / 2 (13) 

where C?(x), the Gegenbauer polynomials, are defined by the serie: 

(l-2xs + s 2 )- a =Y J C?{x)s j (14) 
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Then, one obtains: 

C'Z(x) = 1 (15) 

C?(x) = 2ax (16) 

C%(x) = 2a(a + l)x 2 -a (17) 

C3 (a;) = + l)(a + 2)x 3 - 2a(a + l)x (18) 

C%{x) = \a{a + l)(a + 2) (a + 3)x 4 - 2a(a + l)(a + 2)x 2 + ^a(a + 1) (19) 

From these relations it is possible to obtain all the F n i for the Iodine and Sodium 
atoms. 

Let us now consider the counting rate R.A-e expected in a detector made by various 
kinds of target atoms - as it is the case of Nal(Tl) - identified by the b index: 

R A _ e=NT £^ee {v 2 ) £ £ N nl F^ (q) p 6(m a - E nl ) (20) 
b=Na,I nl 

with q ~ p ~ y // 2m e (m a — E n i). There the second sum is running over nl levels - each 
of them with N n i electrons - allowed by the Heaviside step function 9(m a — E n {). The 
binding energy E n i of Na or I atoms have already been reported in Tab. 1. 

3.1.3 The Primakoff effect 

As regards the Primakoff effect (see Appendix): 

(21) 



where 9 is the scattering angle between a and 7, k is the a momentum and $ is the 
electric potential of the target nuclei. 

Since, the Nal crystal is a ionic crystal, the electric potential generated by each 
ion is well approximated by two components, a Yukawa one shielded by the bound 
electrons and a Coulomb one at long range, respectively: 

* (Zb±l)e -J*. e 

where the upper(lower) sign is for b = I(b = Na) and r& is the screening radius of the 
atom. The Fourier transform of the electric potential is: 



Hence, the counting rate Rp r i m can be expressed as: 
6ir 



R Pnm = NTP f^ Y, IWlV) ( 24 ) 



b=Na,I 



with q ~ m a . 

Therefore, also in this for the axioelectric one, a non-zero modulation term 

is present. 
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3.1.4 Conclusions 

In conclusion, all the kinds of interaction discussed above for the pseudoscalar candi- 
date contribute to the costant part of the signal, while the Compton-like interaction 
does not contribute to the modulation part of the signal. Moreover, as it can be 
easily demonstrated, for the pseudoscalar case the axioelectric contribution to the to- 
tal expected counting rate is largely dominant with the respect to the Primakoff and 
Compton-like on nuclei contributions at least in all the "natural" cases, where gage/me 
is not lower than a factor <~ 10~ 3 the coupling constant to mass ratios of the other 
charged fermions; in addition, it still remains at least one order of magnitude larger 
than the one due to the Compton-like effect on electrons, for a particle mass below ~ 
6 keV. 

3.2 The data analysis and results for the DM pseudoscalar can- 
didate 

A complete data analysis of the 107731 kg • day exposure from DAMA/Nal in this 
framework would offer an allowed multi-dimensional volume in the space defined by 
m a and by all the g a jf coupling constants to charged fermions. Since a graphic rep- 
resentation would be practically impossible, we will present some particular slices of 
such an allowed volume. 

In the following calculations, we have considered some of the possible uncertainties 
in the parameters needed for the counting rate evaluations. In particular, the screen- 
ing radia of Sodium and Iodine atoms have been taken in the range: = (2.0 ±0.5) 
keV and ^_ = (5.0 ± 0.5) keV, respectively. Moreover several halo models have been 
considered [14, 15] either spherically symmetric matter density with isotropic velocity 
dispersion or spherically symmetric matter density with non-isotropic velocity disper- 
sion or axisymmetric models or triaxial models. The parameters of each halo model 
have been chosen taking into account the available observational data; in particular 
the local velocity has been varied within its allowed range: v$ = (220 ± 50) km s _1 
(90% C.L.) and local density phaio has been varied within the range evaluated as in ref. 
[13] taking into account the following physical constraints: i) the amount of flatness of 
the rotational curve of our Galaxy, considering conservatively 0.8 • vq < < 1.2 • vq, 
where v^t is the value of rotational curve at distance of 100 kpc from the galactic 
center; ii) the maximal non dark halo components in the Galaxy, considering con- 
servatively 1 • 10 10 M Q < M vis < 6 • 10 10 M Q [40, 41]. Although a large number of 
self-consistent galactic halo models, in which the variation of the velocity distribution 
function is originated from the change of the halo density profile or of the potential, 
have been considered, still many other possibilities exist. The parameter (Vg) has been 
evaluated for each considered halo model and for the given local velocity value. 

The results presented in the following by using the DAMA/Nal annual modulation 
data over the seven annual cycles are calculated by taking into account the time and 
energy behaviours of the single-hit experimental data through the standard maximum 
likelihood method 5 . 



5 Shortly, the likelihood function is: L = Ilj^e Mi J fc * jfc , , where iVj 7 fc is the number of events 
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In particular, the likelihood function requires the agreement: i) of the expecta- 
tions for the modulation part of the signal with the measured modulated behaviour 
for each detector and for each energy bin; ii) of the expectations for the unmodulated 
component of the signal with the respect to the measured differential energy distri- 
bution. In the following for simplicity, the results of these corollary quests for the 
pseudoscalar candidate particles are presented in terms of allowed regions obtained as 
superposition of the configurations corresponding to likelihood function values distant 
more than 3er from the null hypothesis (absence of modulation) in each of the several 
(but still a limited number) of frameworks. Obviously, these results are not exhaus- 
tive of the many scenarios (still possible at present level of knowledge) for these and 
for other classes of candidates, such as the WIMPs we already deeply investigated 
[6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. 

First of all, as already mentioned, the axioelectric contribution is dominant with 
the respect to the Compton-like and Primakoff effects in all the "natural" cases; thus, 
the results can be presented in terms of only two variables g age and m a . The allowed 
region in the plane defined by these two variables has been calculated considering the 
DAMA/Nal results on the model independent annual modulation signature and has 
been reported in Fig. 3. It is worth to note that, for the reasons reported above, this 
allowed region is almost independent on the adopted g aau and g add coupling constants. 
The allowed region reported in Fig. 3 can only marginally be affected by the results 

-5 

g 10 f 1 




0.5 1 2 3 4 5 6 7 8 910 

m (keV) 

Figure 3: Pseudoscalar case: region in the plane g a se vs m a allowed by the DAMA/Nal 
annual modulation data in the considered model framework. See text. 



collected in the i-th time interval, by the j'-th detector and in the fc-th energy bin. N%jk follows a 
Poissonian distribution with expectation value /lijfc = [bjk + So.fc + ^m.fc ' cosu)(ti — to)]MjAtiAEej^. 
The unmodulated and modulated parts of the signal, Sq ^ and S m k cosuj(ti — to), respectively, are 
functions of the coupling constants and of the light boson mass; bj^ is the background contribution; 
Ati is the detector running time during the i-th time interval; ej^ is the overall efficiency and Mj is 
the detector mass. 
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already presented at low energy by low-background ionization detectors; in fact, due 
to their energy resolution and to their quoted counting rate at low energy, their results 
do not rule out a particles with m a < 3 keV and, for m a > 3 keV, a particles with 
g aSe < 2 x 10- 10 . 

Some strongly model dependent astrophysical limits on the g aSe can be found in 
literature (see e.g. [42]) by studying the globular cluster stars; however, these con- 
straints only apply to particles with masses much below few keV, which is the typical 
core temperature of the stars. 

DAMA/Nal allowed region depending 
on the values of all the g oT , 




m Q (keV) 

Figure 4: DAMA/Nal allowed region at 3cr C.L. in the plane g a ~ n vs m a for a light 
pseudoscalar candidate (crossed hatched region). All the configurations in this region 
can be allowed depending on the values of all the g a ff- Examples of two of the many 
possible models: i) upper black region: coupling only to electrons; ii) lower black region: 
coupling (proportional to the m a ) through the (weak) isospin to quarks and leptons. 
For details see text. The indicative region of the Kaluza-Klein (K.K.) pseudoscalar 
axion credited in ref. [30] from the analysis of indirect observations and the region of 
the DFSZ and KSVZ models (P.Q. axion) are shown as well. The solid line corresponds 
to a particle with lifetime equal to the age of the Universe; at least all the .ga 77 's below 
this line are of cosmological interest. See text. Thus, a pseudoscalar DM candidate 
could also account for the DAMA/Nal model independent result as well as the WIMP 
solutions already discussed elsewhere. 

Let us now investigate the cosmological interest of the allowed a particle, actually 
the case of an a lifetime of order or larger than the age of the Universe. As reported 
in eq. (7), this lifetime is directly connected to the g a ~ n value. Hence, we show - for 
a general case - in Fig. 4 the DAMA/Nal region allowed in the "usual" plane g ail vs 
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m a (crossed hatched region). The upper bound on g 077 is given when the Primakoff 
effect is largely dominant (that is, if g a se = g a uu = g a dd = and only contributions 
from other charged fermions are present). All the other values of ga-y-y below this 
upper bound are allowed depending on the values of all the g a jf. Just for example, 
two particular (of the many possible) sets of g a ^ are reported in Fig. 4 as upper black 
region and lower black region, respectively: the case of the coupling only to electrons 
t g "" a ) and the case of g aUu and g a dd according to the scenario of ref. [26] for 



(ffa77 — 

the Maioron. This last model assumes that q n r f oc rat: = = —S^l- sam c 
relations also hold for the other families. Replacing these assumptions in the eq. (8), 
we have at the first order: 



+ 



if add , 4.9 auu 

3m d 3m u 



, 1 4 
1+ 3-3 







while at the second order: 



127rm P 



1.4 10" 6 GeV' 1 



mi 



■ 9ae 



1 keV 



(25) 



(26) 



Fig. 4 also shows the result of Kaluza-Klein axions credited in ref. [30] from 
analysis of indirect observations; it is well embedded in the DAMA/Nal allowed region. 
Moreover, when considering the existing uncertainties (qualitatively reported in figure 
as shaded area), the result of ref. [30] can also be in agreement with the DAMA/Nal 
region allowed for the scenario of ref. [26] depicted in the lower black region. Moreover, 
the DAMA/Nal allowed region is also comparable with the predictions of ref. [31] of 
axion-like particles as UHECR. 

The solid line in Fig. 4 corresponds to a particle with lifetime equal to the age 
of the Universe; therefore, at least all the configurations with g ail below this line 
are of cosmological interest. Moreover, in principle, it might be possible that the 
configurations above this line would also become of interest in case of some exotic 
mechanism of the a particle production. 

Thus, a pseudoscalar DM candidate could also account for the DAMA/Nal model 
independent result as well as the WIMP solution for several (but still few with the 
respect to the possibilities) corollary model dependent quests for the candidate particle 
already discussed elsewhere [6, 7, 8, 9, 10, 11, 12, 13, 14, 15]. 



4 Scalar case 

The effective Lagrangian of interaction of the scalar particle (we remind that here and 
hereafter the scalar particle is named h) is given by the sum of the Yukawa couplings 
to the fermions and the effective coupling to photons through a charged fermion loop: 

L int = g hf - f h ff + 9 -^h F^F"" (27) 

The corresponding diagrams are similar to those reported in Fig. 2, just replacing 
a with h. Also in this case the scalar h particle can decay into photons with lifetime 
T h = g/ 4 ^ i as regards the decay into pair of fermions similar considerations of sec. 
3 still hold. 
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The effective coupling constant to two photons gu-yy is for run «rrtf [35]: 



9h~yy 



E 



9hJ}Q)<* 
7r rrif 



3 180 mj \mj 



(28) 



4.1 Direct detection of the scalar particle 
4.1.1 The Compton-like effect 

The differential cross section of the Compton-like effect (see also Fig.l) for h particle 
impinging on point-like free fermion is given by the sum of two contributions (see the 
Appendix) : 



da C - 



■J 



8ir mj\v\ 



" l h 



(29) 



The first term is similar than that obtained for the pseudoscalar interaction and, since 
it is proportional to the inverse of h velocity, it does not contribute to the modulation 
amplitude of the signal; the second term is instead proportional to h velocity and it 
does contribute to the modulation part of the signal. 

The expected counting rate, Rc-u hi a detector can be written as the sum of 
Rc-i,eiect and Rc~i,nuci given by the conversion on the atomic electrons and by the 
conversion on the nuclei, respectively. 

For the conversion on electrons: 



Rc-l,e 



N 7 



■ E 

A=Na.I 



Pa aZ A9hee m h 

2mi 



4 mi , , 
3 m 2 h 



(30) 



For the nucleus: 



Rc- 



l,nucl 



= N 7 



E 

v=Na,I 



A,eff9 h AA mh 

2m\ 



4 m A 
3m 2 h 



F 2 



A,Nucl 



(31) 



The scaling law for the g^AA nas been evaluated considering the coherent contri- 
bution of all the quarks: g hAA = A [{g h u U + ^9hdd) + f i.9huu - 9hdd)] ■ Since the ratio 
Z/A is quite similar for the nuclei (0.48 and 0.42 for the 23 Na and for the m I, re- 
spectively), we can write without losing generality g hA A = ^9hNNy where g^NN ^ s the 
effective coupling constant to a single nucleon. The effective charge of nuclei, ZA, e ff, 
takes into account the screening of the inner electrons; a reasonable approximation for 
uih ~ keV is ZNa,eff = 9 and Zi. e ff — 35 (see Tab. 1). The form factors for the 
nuclei can be considered as one since the nuclei are practically point-like for the energy 
of interest. 

Finally, the Compton-like effect provides in the present case a counting rate varying 
along the year. 
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4.1.2 The photoelectric-like effect 

The cross section of the photoelectric-like effect of a scalar h particle can be evaluated 
as reported in the Appendix and is: 

V p -e,nl=9lJ-^\ihF 2 n M) (32) 
|fc| 

The counting rate R p - e can be calculated straighforward as done in sec. 3.1.2 for the 
pseudoscalar case. 

Rp _ e=Nr E^^l £ J2 N ^ F nM)P^(rn h -E nl ) (33) 

h b=Na,I nl 

It is worth to note that the photoelectric-like effect for the scalar case gives contribution 
only to So part of the signal and does not provide any significant time dependent 
contribution. 



4.1.3 The Primakoff effect 

The cross section and the expected counting rate for the Primakoff effect are equal 
than those reported for the pseudoscalar case, with the obvious substitution of the 
coupling constant (see Appendix). As in that case, a non-zero modulation term is 
expected. 



4.1.4 Conclusions 

In conclusion, all the kinds of interactions discussed above for the scalar candidate 
contribute to the costant part of the signal, while the photoelectric- like interaction does 
not contribute to the modulation part of the signal. However, the photoelectric-like 
interaction gives a dominant contribution with the respect to the Compton-like process 
on the electrons and to the other effects. Therefore, the coupling to the electrons does 
not produce any significant time variation of the signal. Thus, in the following we will 
investigate the case of a scalar h particle coupled only to the hadronic matter; in this 
case, non-zero modulation term of the signal is expected by the contributions of the 
Compton-like effect on nuclei and of the Primakoff effect. 



4.2 The data analysis and results for the DM scalar candidate 

Processes involving DM scalar candidate coupled to electrons do not contribute to the 
annual modulation of the signal. Moreover, the upper limit on ghee, which can be eas- 
ily derived from the energy distribution measured by DAMA/Nal (given elsewhere), 
ranges from ~3x 1CP 16 to ~ 10~ 14 for mu between ~ 0.5 and 10 keV. Thus, in the 
following we investigate the case of a scalar h particle coupled only to the hadronic 
matter (ghee « 9huu,ghdd)- This framework would offer an allowed multi-dimensional 
volume in the space defined by m/, and by all the g h jf coupling constants to quarks. 
However, to account for cosmological interesting h particle lifetimes, only the con- 
figurations with negligible Primakoff effect contribution can be investigated and only 
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the diagrams at the tree level with quarks are considered. Therefore, in this case the 
Compton-likc effect on nuclei is the major contribution to the cross section. Hence, 
according to the scenarios already described in sec. 3.2, we report in Fig. 5 the region 
allowed by the DAMA/Nal data in the plane g^NN vs m h- This allowed region can 
only marginally be affected by the results already presented at low energy by low- 
background ionization detectors, which do not rule out h particles with < 3 keV 
and, for m/j > 3 keV, h particles with g^jvAf < 8 x 1CP 7 . Finally, the strongly model 
dependent globular cluster constraints [43] only apply to particles with masses much 
below few keV. 




0.5 1 2 3 1 5 6 7 8 910 



m h (keV) 

Figure 5: DAMA/Nal allowed region at 3cr C.L. in the plane g^NN vs mh for a scalar 
light boson DM candidate coupled only to the hadronic matter in the given frameworks. 
To account for cosmological interesting lifetimes, here only the configurations with 
negligible Primakoff effect contribution have been selected. See text. 

Obviously, the region in Fig. 5 does not allow a direct information about the effec- 
tive coupling constant to photons, 37177, and, therefore, about h lifetime; however, the 
large number of free coupling constants allows to expect the existence of a large num- 
ber of h configurations of cosmological interest. For example, assuming as reasonable 
hypothesis that all the ghqq's would be of the same order of magnitude, Th would be 
dominated by the quarks of the first family (see eq. (28)). In this particular case, we 
can report in the g hau vs g h g d plane the configurations of cosmological interest allowed 
by the results of DAMA/Nal. Fig. 6 shows the shaded regions obtained by requiring 
that the configurations: i) arc allowed in the plane g^NN vs m h a $ reported in Fig. 
5; ii) TO/i > 0.3 keV; iii) are of a cosmological interest (t^ > the age of the Universe). 
The uncertainties on the first family quark masses have also been accounted for [43]: 
m u ~ 3.0 ± 1.5 MeV, m d ~ 6.5 ± 2.0 MeV. 

It is worth to note that we have not included in the presented calculation a possible 
contribution due to the scalar interaction of h particle to the gluons, which can natu- 
rally explain the pure coupling to the nuclear matter. In principle, this contribution 
might increase the configurations of cosmological interest and allowed by DAMA/Nal. 

Finally, also a scalar DM candidate could account for the DAMA/Nal model in- 



1G 



Figure 6: DAMA/Nal allowed configurations of cosmological interest in the ghuu vs 
9hdd plane obtained by requiring that they: i) are allowed in the plane g^NN vs m h 
in Fig. 5; ii) rrih > 0.3 keV; iii) Th > the age of the Universe. The uncertainties 
on the first family quark masses have also been accounted for: m u ~ 3.0 ± 1.5 MeV, 
m,d ~ 6.5 ± 2.0 MeV. In the inset: a magnification of the regions around (0,0). See 
text. Thus, also a scalar DM candidate could account for the DAMA/Nal model 
independent result as the previous case of pscudoscalar DM candidate and the WIMP 
solutions already discussed elsewhere. 

dependent result as the previous case of pseudoscalar DM candidate and the WIMP 
solutions already discussed elsewhere. 



5 Conclusions 

In this paper we have discussed in some details another class of DM candidate par- 
ticles: the pseudoscalar and scalar light bosonic candidates; particular care has been 
devoted to the study of the processes for their detection in a suitable underground ex- 
perimental set-up. In future we also plan to explore additional scenarios which can - 
in principle - exist such as a particle with mixed pseudoscalar and scalar interaction, a 
triplet of (pseudo-) scalar particles (which can "naturally" allow a quasi-stable bosonic 
candidate), particles with gluon coupling, spin 1/2 light dark matter candidates such 
as those of ref. [49, 50], etc.. 

Moreover, it is worth to explore also the possibility that (electron coupled) pseu- 
doscalar a particles produced by high energy sources in the galactic bulge would - by 
interaction cither with interstellar medium or with relic a particle - also contribute to 
positron production and to the 511 keV annihilation 7 rays from the bulge of Galaxy 
observed in ref. [44]. Finally, just candidates produced in the early Universe and 
decoupled from the primordial Universe have been considered here; other mechanisms 
of light bosons production can in principle be considered, such as their production by 
a heavier particles decays or "exotic" cosmological evolution of the Universe, etc.. 
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In conclusion this paper, far to be exhaustive on the topics (e.g. many other 
scenarios can be considered for the galactic halo), has however shown that - in addition 
to the WIMP cases already discussed by DAMA collaboration elsewhere ([14, 15] 
and references therein) - there is also possibility for a bosonic candidate with axion- 
like phenomenology and with mass < 6-7 kcV to account for the 6.3 a C.L. model 
independent evidence for the presence of a particle DM component in the galactic halo 
observed by DAMA/Nal. 

The new higher sensitive DAMA/LIBRA set-up now in operation deep under- 
ground at the Gran Sasso National Laboratory of I. N.F.N, will allow to restrict the 
possibility in the corollary quests for the DM candidate particle as soon as it will have 
collected a well competing exposure with the respect to the previous DAMA/Nal. 



6 Appendix 

6.1 Compton-like cross section 
6.1.1 Pseudoscalar case 

Let us consider the diagram of Fig. 1 where a pseudoscalar a particle with quadrim- 
pulse k is absorbed by a charged fermion of quadrimpulsc p with the subsequent 
emission of a photon with quadrimpulse k' . The outgoing fermion has quadrimpulsc 
p' . The matrix element for Compton-like conversion is given by the sum of the contri- 
butions of this diagram and of the "crossed" one. We get: 



M = ieQ f g af - f U(p') 



ft + mf 



(k+pf 



m. 



'75 + 75 



ft + m f 



i f (p - k 1 ) 2 - mj 

Using the standard algebra of the 7 matrices, the square of matrix element averaged 
on the initial fermion polarizations and summed on the final electron polarizations can 
be written as: 



U(p) 



(34) 



| M p = ^ZlLTr [(/ + m f ) (^fa - {j + m } ) - g£ 



(35) 



After some algebric manipulations and summing on the outgoing photon polarizations, 
we obtain: 



|M| 2 = e 2 Qjg 2 aff 2A + B\k\ 2 sin 2 {e) 
where: 

A = 

and 



2k' -p 



2p-k 



- in 



2p-k 



B = -- 



+ 



2p-k' 
8(k-p) 



a _ 2 



m 2 , +2p ■ k (m 2 a + 2p ■ k) 2 
The cross section is given by (see for example [45]) 

d 3 k' 



da = 



\M[ 



1 1 



d 3 p' 



\vi - v 2 \ 2p 2k (2tt) 3 2 k' (2tt) 3 2 p' 



7 (2tt) 4 <5V + k' -p-k) 



(36) 



(37) 



(38) 



(39) 
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For mf » m a and considering that both the a particle and the fermion can be treated 
in the non-relativistic limit {E 1 ~ E a ~ m a ), we obtain A — > ^f- and _B|fc| 2 — » — -^j- ~ 

2 / / 

— /3 2 (/? — 10 -3 is the a particle velocity in the galactic halo). Therefore, the cross 

section can be evaluated by considering that the A term is dominant in the squared 
matrix element: 



da 



ff m a 



dtt dE^ 



TUfTT 



l f 8\k\ 



S(E 1 - E a ) 



(40) 



and eq. (9) can be inferred. 



6.1.2 Scalar case 



The matrix element for Compton-like conversion for scalar case can be written, using 
the same formulation as in the previous case, as: 



M = eQ f g hff U(p') 



U(p) 



(41) 



~(k+p) 2 — m 2 - {p—k') 2 —mf 

and the square of matrix element averaged on the initial fermion polarizations and 
summed on the final electron polarizations: 

W = ^¥^Tr [(/ + m f ) 4 - {i> + mf) - ^) i\ (42) 

After some algebric reductions and summing on the outgoing photon polarizations, we 



obtain: 



\M\ 2 
where: C = B + 



16ml 



(m£+2p-fc) : 



2 A + C\k\ 2 sin 2 (9) (43) 
-. For m,f » m a and in the non-relativistic limit: C\k\ 2 — ► 



= Aj3 2 . Hence, on the contrary of the pseudoscalar case, both the A and C 
terms in the squared matrix element can contribute at different extent to the total 
cross sections depending on the h particle mass and on the considered fermion, either 
electrons or nuclcons. Therefore, the non-relativistic cross section is: 



da 



aQ)g 



2 „2 _ 
hff 



dE^dn 



87T/3TO 2 - 



m. 



5{E 1 - E h ) 



(44) 



from which eq. (29) can be derived. 



6.2 Photoelectric-like cross section 
6.2.1 Pseudoscalar case, axioelectric process 

From the Lagrangian (6) it is possible to obtain in the non-relativistic limit the 
Schroedinger-Pauli equation for the electron in presence of the a field and of the 
electromagnetic field: 

(v + iel) a -Via) 

(id + eA ) V + ^— '—^ - —a ■ (V x A)ip - g a - ee J ( > $ = (45) 



2m P 



2m e 



2m e 
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It is evident, therefore, as in the non-relativistic limit, the a particle interaction to 
fermions is spin-dependent. 

The cross section of the process can be evaluated through the Fermi golden rule: 



da 



Ji 



2tt^ 



\M, 



if\ 



J, 



-Pf 



(46) 



when considering the perturbation term in the electron Hamiltonian (see for example 

[46]),K=5ae-e^. 

Considering the a field and the wavefunction of the escaping electron as plane 
waves and the wavefunction of initial bounded electron in the nlm shell as ip n im, the 
transition matrix element can be written as: 



\M, 



-ip-x 



-gaeeCloe 



2 m, 



(47) 



The incident flux is: Ji — 
outgoing electron is: 
Vp 2 dp 



-i(a*Va — aVfl*) = 2a^k and the density of states for the 



Pf 



3 dCl6(E a - E nl - E e ) = me J P \ d 3 Ee dnd(E a - E nl - E e ) (48) 



(2tt) 3 — -^-» (2tt) 3 

where E n i is the binding energy of the atomic electron in the nl shell. Summing over 
the final electron polarizations and averaging on the initial polarizations, we get: 



da 



dE e dn 



= 9le 



k 



4m P 



(27T) 



-5(E a — E n i — E e ) 



(49) 



From this formula, eq. (12) can be deduced straighforward. 



6.2.2 Scalar case 

The procedure to evaluate the cross section for the scalar photoelectric-like case is 
similar to that applied before for the pseudoscalar case. Considering the Lagrangian 
given in eq. (27), the last term of the Schroedinger-Pauli equation (45) has to be 
replaced by the term gheehip = —Vuip. Therefore, the cross section for the process in 
the non relativistic limit can be written as: 

2 



da 



= 9 2 hi 



I'- 



i(p-k)x 



dE e dVL ™ ee k 
and the eq. (32) can be inferred. 



P 



(27T): 



;S(Eh - E n i - E e ) 



(50) 



6.3 Primakoff cross section 
6.3.1 Pseudoscalar case 

Let us now consider the case of the Primakoff conversion of the pseudoscalar a particle 
of momentum k into a photon of momentum p in the static atomic field F pa (x). The 
S-matrix element is: 

2 



\S 



d 4 x^ae^F^F pa (x) 



(51) 
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where the factor 2 takes into account the two photon crossed diagrams. Let us now 
consider the a field and the wavefunction of the outgoing photon as plane waves: a = 



a e 



-ik-x an( } a _ f 



The cross section is given by (see for example [47]): 



For the static electric field: e^ pa F prJ (x) = 2e^ M E l (x) = 



da = V 



d 3 P \S ft f 
(2irf T Ji 



(52) 



-i(a*\7a — aVa*) 



2a\k 



As for the previous case, the incident flux is given by: J. L 
and the time integral for a transition in the time T can be obtained following the 
prescriptions given in rcf. [47] . Summing on the possible polarization of the outgoing 
photon, we obtain the differential cross section of the process: 



da 



Q 2 

dVLdE 1 ~ 167?' 



2 (e)E* k 



d 3 x^(x)e^ s 



S(E 7 - E a ) 



(53) 



Similar results have been obtained in ref. [48]. It is worth to note that this case 
also includes the ultrarelativistic limit, considered for example for the solar axion 
investigations. Eq. (21) can be directly achieved by the previous eq. (53). 



6.3.2 Scalar case 

For a scalar h particle the S-matrix element for Primakoff conversion is: 



\Sfi\ 



(54) 



Using the same procedure as for the pseudoscalar case and reminding that the electric 
static field is: 

F^{x) = - [r^rF - (55) 
one can obtain the differential cross section as for the pseudoscalar case: see cq. (53). 
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